Josephson current and 7r— state in ferromagnet 
with embedded superconducting nanoparticles 

A. V. Samokhvalov (1) , A. I. Buzdin (2) 
W Institute for Physics of Micro structures, 
Russian Academy of Sciences, 603950 Nizhny Novgorod, GSP-105, Russia 
( 2 ) Institut Universitaire de France and Universite Bordeaux I, 
CPMOH, UMR 5798, 33405 Talence, France 

Abstract 

On the basis of Usadel equations we investigate superconductor/ferromagnet/superconductor 
(S/F/S) hybrid systems which consist of superconducting nanostructures (spheres, rods) embedded 
in ferromagnetic metal. The oscillations of the critical current of the S/F/S Josephson junctions 
with the thickness of ferromagnetic spacer between superconducting electrodes are studied. We 
demonstrate that the n state can be realized in such structures despite of a dispersion of the 
distances between different parts of the electrodes. The transitions between and ir states at some 
thickness of ferromagnetic spacer can be triggered by temperature variation. 

PACS numbers: 
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I. INTRODUCTION 



The particularity of the proximity effect in superconductor/ferromagnet (S/F) hybrid 
structures is the damped oscillatory behavior of the Cooper pair wave function inside the 
ferromagnetic (for the reviews see^). In some sort it is a manifestation of the Larkin- 
Ovchinnikov-Fulde-Ferrell (LOFF) state induced in ferromagnet (F) near the interface with 
superconductor (S). In contrast with original LOFF, which is possible only in the clean 
superconductors, the damped oscillatory S/F proximity effect is very robust and exists also 
in the diffusive limit. This special type of the proximity effect is at the origin of the ir 
Josephson S/F/S junction^ which has at the ground state the opposite sign of the super- 
conducting order parameter in the banks. Firstly the ir- junction was observed at experiment 
in Ref.- and since then a lot of progress has been obtained in the physics of ir- junctions and 
now they are proven to be promising elements of superconducting classical and quantum 
circuits^. Different manifestations of unusual proximity effect and the ir states have been 
observed experimentally in various layered S/F hybrids^-. The proximity induced switching 
between the superconducting states with different vorticities in multiply connected hybrid 
S/F structures was suggested recently in^*ii. Theoretical studies and experiment sf£~— both 
demonstrated that in the diffusive limit the spin-flip and spin-orbit scattering lead to the 
decrease of the decay length and the increase of the oscillating period. In addition the spin- 
flip scattering may generate the temperature induced transition from to ir state of the 
junction^ 1 ^ . 

Naturally at the first stage of the work on the S/F/S junctions the systems with planar 
geometry and well controlled F-layer thickness were considered. However, now when the ir 
state proven to be very robust vs different types of the impurities scattering^^ (magnetic 
and non-magnetic), interface transparency^ 1 ^ it may be of interest to address a question 
how the ir junction could realize in the S/F/S systems with a bad defined thickness of F- 
spacer, in particular for two superconducting particles imbedded in ferromagnet or between 
a flat superconducting electrode and a small superconducting nanoparticle (such situation 
could be of interest for the STM-like experiments with a superconducting tip). This question 
is non-trivial because the transition from to ir state occurs at the very small characteristic 
length^ = y/Df/h, where Df is the diffusion constant in ferromagnetic metal and h is the 
ferromagnetic exchange field, and the typical values of are in the nanoscopic range. We 
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could expect that when the variation of the distance between different parts of S-electrodes 
is of the order of £/ the ir state would disappear. Our calculations show that it is not the 
case and once again the tt state occurs to be very robust and the transition between and 7r 
states is always present at some distance and also can be triggered by temperature variation. 

In this paper we present the results of a theoretical study of the peculiarities of the 
proximity effect and Josephson current in S/F hybrids which consist of superconducting 
nanostructures placed in electrical contact with a ferromagnetic metal. The paper is orga- 
nized as follows. In Sec. II we briefly discuss the basic equations. In Sec. Ill we calculate 
the Josephson current in two model hybrid S/F/S systems. The first system consists of two 
superconducting rod-shaped electrodes imbedded in ferromagnet. The second one is a S/F 
bilayer with a superconducting spherical particle at the surface of the ferromagnetic layer. 
We examine the temperature dependence of the critical current of S/F/S junction between 
the flat superconducting electrode and the S-particle taking into account the spin-flip scatter- 
ing. For both cases the S/F interface transparency between superconducting nanoparticles 
and ferromagnet is assumed to be low to prevent from superconductivity destruction due to 
proximity. We summarize our results in Sec. IV 

II. MODEL AND BASIC EQUATIONS 

Since the models of S/F/S junctions we are going to study consist of superconducting 
particles embedded in a ferromagnetic matrix or placed on a ferromagnetic substrate, we 
start from a description of the damped oscillatory behavior of the Cooper wave function 
induced by such particle in a ferromagnet. 

We assume the elastic electron-scattering time r to be rather small, so that the critical 
temperature T c and exchange field h satisfy the dirty-limit conditions T c r 1 and hr <C 1. 
In this case a most natural approach to calculate T c is based on the Usadel equations^ for 
the averaged anomalous Green's functions Ff and F s for the F- and S-regions, respectively 
(see^ for details). These equations are nonlinear but can be simplified when the tempera- 
ture is close to the critical temperature T c or at any temperature of the F-layer when the 
transparency of S/F interface is low. In the F- region the linearized Usadel equations take 
the form 

- ^ V 2 F f + (\cj\+ihsgncu + 1/t s ) F f = , (1) 
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where Df is the diffusion coefficient in ferromagnetic metal and u are the Matsubara fre- 
quencies, to = 2irT(n + 1/2), and r s is the magnetic scattering time. We consider the 
ferromagnet with strong uniaxial anisotropy, in which case the magnetic scattering does not 
couple the spin up and spin down electron populations. Restricting ourselves to the case of 
superconducting inclusions with cylindrical or spherical symmetry (cylindrical rod-shaped 
or spherical particles of radius R s ), one can easily find the following solutions of Eq. (JTJ), 
which describe the distribution of anomalous Green's function Fj in ferromagnet (r > R s ) 
surrounding a superconducting cylinder 

F c f (r) = AK (qr) , (2) 

or superconducting sphere 

Fj(r) = Aexp(—qr)/qr, (3) 
where K (z) is the Macdonald function, 



q = ^2/D f ^\uo\+ihsgnu + l/r s (4) 

is the characteristic wave number of the order parameter variation in the F-metal, and the 
amplitude A is determined by the boundary conditions at the S/F interface (r = R s )r^: 

a s d r F s = a n d r F f , F s = F f - ^ n d r F f . (5) 

The S/F interface between a particle and ferromagnet is assumed to be characterized by 
the dimensionless parameter 7^ = Rb(J n /^n related to the boundary resistance per unit area 
Rb. Here = -J ' D s ^/2ttT c is the superconducting (normal-metal) coherence length, 
a s and a n are the normal-state conductivities of the S- and F-metals, and d r denotes a 
derivative taken in the radial direction. We will assume that the rigid boundary condition 
7b 3> min{£ s <T n /£ n er s , 1} is satisfied, when the inverse proximity effect and the suppression 
of superconductivity is S-metal can be neglected^ 1 ^. As a result, the pair amplitude F s (r) 
at the S/F interface is equal to the one far from the boundary: 

F S (R S ) = ^L= = ^-G n , (6) 

where A is the superconducting order parameter, and 



is normal Green's function. Using the solution ([2} or (j3J) we obtain from the Eqs. (jSJ) and ([6]) 
the value of Ff at the S/F boundary r = R s , and the amplitudes A for both cases. Finally, 
the expressions 



n J ftr (g/? s )+ 7fe en.g^i(^. 



W = - G- ., ■ , 71 ■ i;n u • (9) 



Z7 n r(l + 7^„.(g+l/i? s )) 
describe the the damped oscillatory behavior of anomalous Green's function Ff in ferromag- 

net surrounding the superconducting cylinder or sphere, respectively. 

The general expression for the supercurrent density is given by 

£ = ^ £ (*> v *> - F f vp f) > ( 10 ) 

u>, cr=± 

where Ff(r,u>) = FUr, —oj). 



III. CRITICAL CURRENT OF JUNCTIONS WITH SUPERCONDUCTING PAR- 
TICLES 

Now we proceed with calculations of the Josephson critical current for two examples of 
mesoscopic hybrid S/F systems. The first one is two identical superconducting cylindrical 
rod-shaped electrodes surrounded by a ferromagnetic metal (see Fig. [I]). The second one is 
a S/F bilayer with a superconducting particle at the surface of the ferromagnetic layer (see 
Fig. ED. 



A. S/F/S junction between two superconducting rod 

Consider two superconducting cylinders of a radius R s embedded in ferromagnet as it 
is shown in Fig. [U The distance between the cylinder axes is d > 2R S . These rod-shaped 
electrodes form a Josephson junction in which the weak link between two superconductors 
is ensured by ferromagnetic neighborhood. The supercurrent 

I s (tp) = I c sin(y2) (11) 

flowing across this structure depends on the phase difference if between the order parameters 
of the rods: 

A li2 = Ae ± ^ /2 . (12) 
5 
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FIG. 1: (Color online) Schematic representation of the F/S hybrid system under consideration: 
two identical superconducting cylindrical rod-shaped electrodes of radius R s surrounded by a fer- 
romagnetic metal. The axes of superconducting cylinders are assumed to be parallel. Figure shows 
the cross section of the structure by the plane (x, y) perpendicular to the cylinder axis. 



For large enough distance between the superconducting cylinders (a = d — 2R S > 2£f), the 
decay of the Cooper pair wave function in ferromagnet in the first approximation occurs 
independently near either of the electrodes and can be described by the solution (JSJ). There- 
fore the anomalous Green function Ff(r) in ferromagnet nearby the plane x = may be 
taken as the superposition of the two decaying functions (jSJ), taking into account the phase 
difference 

j? i , K (qr + )e^ 2 + K (qr^e-^ 2 

f[X,V) ~ u * Ko(qRs)+l b ^qK 1 (qR s ) ' ^ 



where r± = \f(x ± d/2) 2 + y 2 . Using the expression (II 01) . we obtain the sinusoidal current- 
phase relation (jlip in the S/F/S Josephson junction between two superconducting rod- 
shaped electrodes for the case of low transparent S /F interfaces. For the critical current of 
such Josephson structure, we have 
2nTa n ^ A 2 



u)>0 



x Re 



a q 



[KoiqR^ + ^qK^qRs)}' 



dy 



r 



(14) 



where r = \/y 2 + d 2 /A. In the limit of large R s 3> £/ a curvature of the electrodes is not 
essential and the formula (fT4j) coincide with the corresponding expressions for the critical 
current of S/F/S layered structures with a large interface transparency parameter % previ- 
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FIG. 2: (Color online) Influence of the electrode radius R s on the dependence of the critical current 
I c (fl5|) on the distance a between two superconducting rod-shaped electrodes embedded in F-metal 
(1/t s = 0). The numbers near the curves denote the corresponding values of the radius R s in the 
units of The inset gives the zoomed part of the I c (o) line , marked by the dashed box. 

ously obtained in Refs.— . The critical current equation ( fT4l) can be simplified for h irT c 
and R s ,tq ^> £/ and may be written as 

-2(y/y*+cP/4-R s )/t f 

(15) 




tanh — . (16) 

Note that our approach is valid for large enough distance between the superconducting 
cylinders a > 2£j and the first thansition into the it state at a$ is described only qualita- 
tively. It has been demonstrated^ that for the planar S/F/S junction with low interface 
transparency the first transition into the tt state occurs at F layer thickness smaller than 
£/ ( its actual value depends on the the exchange field h and transparency parameter jf, 
).The similar situation is expected for the embedded superconducting particles and to find 
the corresponding interparticle distance we need to solve our problem exactly. 

The dependence of the critical current I c as a function of the distance a between the 



superconducting cylindrical electrodes calculated from Eq. (fT5|) is presented in Fig. [2] for 
several values of the radius R s . From the figure, we see that with increasing the distance 
a, the S/F/S junction undergoes the sequence of O-ir and n-0 transitions when the value of 
I c changes its sign from positive to negative and vice versa. We may roughly estimate that 
the first transition from to it state in S/F/S junction formed by rod-shaped electrodes 
occurs at the thickness of F-layer ao ~ 0.5£/, similar to a S/F/S junction with a low S/F 
interface transparency in the ordinary layered geometry^^. We observe that the distances 
a corresponding to 0-tt and ir-0 transitions grow slightly with decrease of the cylinders radius 
R s due to a dispersion of the distances between different parts of the electrodes (see the inset 



B. S/F/S junction in S/F bilayer with a superconducting particle 

As a second example we consider ferromagnetic film of a thickness dona superconducting 
plate with a transparent S/F interface. The S/F/S Josephson junction is assumed to be 
formed between the flat superconducting electrode and a small superconducting half-sphere 
of radius R s embedded into ferromagnet, as it is shown in Fig. [3J The center of the sphere 
is placed at the surface of the F-film. As before, there is the tunnel barrier (7^ 3> 1) at the 
S/F interface between the superconducting particle and ferromagnetic metal. 

Since we consider a S/F bilayer with a transparent interface then the complete nonlinear 
Usadel equation in the F-layer has to be employed. Using the usual parametrization of the 
normal and anomalous Greens functions Gf = cosGj and Ff = sinO/, the Usadel equation 
is written as 



Note that Eq. ( ITT]) transforms into the linear equation ([TJ in the limit of small 8/ < 1. 
For simplicity we restrict ourselves to the case of thick F-layer (d ^> £/) then the decay 
of superconducting order parameter occurs independently near each S/F interface. In that 
case, the behavior of the anomalous Green's function near each interface can be treated 
separately, assuming that the F-layer thickness is infinite. Following Ref.-^^, the analytical 
solution of the equation ( ITTj) for flat transparent interface at z — d can be written as 



in Fig. [2]). 




(17) 



y/l-e 2 sm 2 (Q f /2) - cos(e//2) 



(18) 



y/1 -5 2 sin 2 (e//2) + cos(G//2) 
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FIG. 3: (Color online) Schematic representation of the S/F/S Josephson junction between the flat 
superconducting electrode and a small superconducting spherical particle of radius R s embedded 
into ferromagnet. Dashed line shows the cross-section of the paraboloid T by the plane (p, z) (|24D : 
p 2 = x 2 + y 2 . Here n is a unit vector along the normal to the surface Y. Figure shows the cross 
section by the plane. 

where 

e 2 = (l/r s ) (\u | + ihsgnu + l/r s ) _1 . 

The integration constant fo should be determined from the boundary condition at the surface 
z = d. As before the rigid boundary conditions is assumed to be valid at z — d: 

Q f (d) = arctan— . (19) 

From Eqs. (JUD,® we get 

(1 " e 2 )F 2 



— , (20) 

(21) 



v/(l-e 2 )F2 + l + l 
IAI 



n 



uj + a/w 2 + |A| 2 

Linearizing the solution ( ITS]) for Gy 1 we obtain the anomalous Green's function in 
ferromagnet (0 < z < d) induced by flat superconductig electrode: 

4 F 

f ~ n — e^"^ . (22) 

v / F^!TT + i 

9 



The total anomalous Green's function in F- layer far from both the S/F interfaces may be 
taken as superposition of the two decaying functions (Q, fl22|) . taking into account the phase 
difference in each superconducting electrode 

A R p-<?( r - R s)+i</V 2 AT? p q(z-d)-itp/2 

f u> n r{l+ lb Uq + ^/Rs)) v/(l-£ 2 )F2 + l + l' 1 ; 

where r = 'x 2 + y 2 + z 2 . 

To derive the general expression for the critical current I c we have to calculate the total 
Josephson current flowing through a virtual surface V: the points of the surface V are 
equidistant from both electrodes of the junction. The surface T is a paraboloid which form 
is described by the equations: 

z = z c -(x 2 + y 2 )/4z c , z c = (R s + d)/2. (24) 

Using the solution (123]) and Eq.f fTOj) . one can arrive at a sinusoidal current-phase relation 
(jTTj) with the critical current (see |A] for details): 

/ J Z ' R * T R C V [ Fn ( A ° n/ U) Q ~ q{d ~ Ra) 

^ [V(l - e 2 )F 2 + 1 + !][!+ lb Uq + 1/R.)] 



n=0 




(25) 



where I = 64ir T c a n £f/e. 

The dependence of the critical current I c fT25l) as a function of the thickness a = d — R s of 
the ferromagnetic spacer separating the superconducting plate and the particle is presented 
in Fig. H]for several values of the particle radius R s and the magnetic scattering time r s . It 
is clearly seen from Fig. H]that with a decrease of the particle radius R s , the position a® of 
the first zero of the critical current is shifted towards larger values of the distance a between 
superconducting electrodes. Figure Hb demonstrates the influence of magnetic scattering on 
the proximity effect and the critical current in S/F bilayer with the particle: decrease of 
the magnetic scattering time r s leads to the decrease of decay length and increase of the 
oscillation period of the anomalous Green's function Ff^. This results in much stronger 
decrease of the critical current in the S/F/S junction with increase of the thickness a, if the 
magnetic scattering time r s becomes relatively small t^ 1 > h. 

Figure [5] shows the temperature dependence of the S/F/S junction critical current I c ( 125]) 
at several values of the thickness of the ferromagnetic spacer between the superconducting 
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FIG. 4: (Color online) The dependence of the critical current I c (|25p on the distance d—R s between 
the superconducting plate and the particle of radius R s embedded in F-metal for different values of 
the radius R s and the magnetic scattering time t s (T/T c = 0.5, h = 3ttT c , 7^ = 10): (a) hr s = 100; 
(b) hr s = 0.5. The numbers near the curves denote the corresponding values of the radius R s in 
the units of £/. 

electrodes. The F-spacer thickness a is chosen close to the first transition from to it state: 
a ~ a — 3£j. The nonmonotonic dependences I C {T) demonstrate - it transition due to a 
change of the temperature. The transition temperature T* (I C (T*) = 0) seems to be very 
sensitive to the size of the superconducting particle. It should be noted however that the 
temperature T* is determined rather by the scale a = d — R s then by the scales d or R s , 
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FIG. 5: (Color online) (a) The dependence of the critical current I c (I25D on the temperature T 
for different values of the radius R s /£,f = 0.95, 1.0, 1.05 and the fixed F-spacer thickness a = 3£/. 
(b) The dependence of the critical current I c (|25p on the temperature T for different values of the 
F-spacer thickness a/£/ = 2.9, 3.0, 3.1 and fixed radius R s = £f. The calculation parameters are 
h = 3ttT c , hr s = 0.5, 76 = 10. 

separately. 
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IV. CONCLUSION 



To sum up, we have analyzed the Josephson effect in S/F/S hybrid structures with a bad 
defined thickness of F-spacer. As an example, we have calculated the Josephson current 
between two rod-shaped superconducting electrodes embedded in ferromagnet or between 
flat superconducting electrode and the small superconducting nanoparticle at the surface of 
the F-layer. For the both cases we have demonstrated the possibility of the realization of tt 
junctions in such hybrid systems. We have studied dependence of the transitions between 
and 7r states both on the size of superconducting particles and the temperature. The ir 
state has been proven to be very robust with respect to a geometry of the S/F/S junction. 
In the dirty limit the transition into n state is determined rather by the thickness of the F- 
spacer between superconducting electrodes then by a shape of the electrodes. Naturally our 
calculations can be easily generalized to the different shape of the S particles (for example 
spherical) with similar conlusions. 

A set of the superconducting particles embedded in a ferromagnetic matrix realize a 
Josephson network. Depending on the geometry of this network and its state (0 or tt) it 
may reveal a spontaneous current similar to that observed in superconducting arrays of tt 
junctions^ . For example the equilibrium phase difference for triangular 2D tt junctions 
network is equal to 2ti/3 which corresponds to the current state. For typical parameters 
Nb/CuNi hybrid system (T c = 9K, £y ps 2nm, p n = l/a n ~ 60/if2cm^ one can get from 
( 1251) the following estimate of the Josephson energy Ej = (PqI c /2ttc of the S/F/S junction: 
Ej/T c ~ 10 4 (/ c //o) < 1, i-e. an observation of spontaneous currents near T c is expected 
to be masked by strong temperature fluctuations. Despite of this restriction we believe 
that intrinsically-frustrated superconducting networks induced by the proximity effect can 
be experimentally observable in such S/F/S composites. In particular, a two-dimensional 
Josephson network of 7rj unctions may serve as a laboratory to study the phase transitions 
with continuous degeneracy^. 

The possibility to fabricate the regular 2D and 3D arrays of Josephson tt junctions and 
monitor the transitions between and tt states simply varying the temperature open inter- 
esting perspectives to study a very reach physics of different phase transition in such systems 
due to interplay between fluctuations, frustration, disorder and dimensionality. 
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Appendix A: Josephson current in S/F bilayer with a superconducting particle 



The general expression for the supercurrent is given by the Eq. fTTUl) , where the anomalous 
Greens function Ff nearby the surface T (see Fig. [3]) may be written as 



F f = A 1 e gz - iip/2 + A 2 e- q ^+** +i<p/2 



A 1 



4F„ e~ qd 



(AG»/«)e* 



(Al) 
(A2) 



y/(l ~ S 2 )F 2 + 1 + 1 l + ^ n (q+l/R s ) 

where r 2 = p 2 + z 2 , and the functions G n , F n are determined by expressions (JTj) and (T2TT) . 
respectively. At the surface T the function Ff and the projection of the vector Vi 7 / along 
the normal 



ii 



Po 



2z c 



to the surface are 



F/| r = ^Ai e qZc - iip/2 + A, 



4z c R s 
p 2 + Az 2 c 



3 qz c +iip/2 \ e -qp 2 /iz c 



(A3) 



(VF f ,n)\ r 



x 



Az r 



2z c e - qp2 / 4z * 

~7p w +V 



-qz c +i<p/2 



(A4) 



> 2 + 4^ c 2 V P 2 + ^4, 
Substitution of Eqs. ( 1A3|) . (1A4|) into the expression for the supercurrent (fTUjl and taking into 
account the symmetry relations q{—oS) = q*(oo), v4 lj2 (— u) = A\ 2 (uj) leads to the following 
formula 



Up) 



327rTa r , 



LU>0 



z 2 R s A 1 A 2 
(p 2 + 4^ c 2 ) 3 /2 



q + 



42, 



Az 2 



J c (p) simp 

-qp 2 /4:Z 



(A5) 
(A6) 
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Further integration of supercurrent density (1A6j) over the surface T 



2z c 



i c = J ' dsj c = j J dppy/p + te* Up) . (A7) 

r o 

results in the expressions (|25p for the critical current I c of the S/F/S Josephson junction 
between superconducting plate and particle. 
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